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Abstract 

The fully self-consistent non-equilibrium Green functions (NEGFs) approach to the quantum 
transport is developed for the investigation of one-dimensional nano-scale devices. Numerical cal- 
culations performed for resonant tunneling diodes (RTDs) of different designs and at different 
temperatures show reasonable results for the potential and electron density profiles, as well as 
for the transmission coefficient and the current-voltage characteristics. The resonant behavior is 
discussed in detail with respect to the quantum-well width, the barrier thickness, and the tempera- 
ture. It is also shown that within the framework of approach used the current noise spectral density 
can be straightforwardly calculated for both the coherent and the sequential tunneling models. In 
qualitative agreement with experiments, obtained results highlight the role of charge interaction 
which causes a fluctuation of density of states in the well and therefore a noise enhancement in the 
negative differential conductance region. 

PACS numbers: 73.63.Hs, 72.70.+m, 73.23.Ad 
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I. INTRODUCTION 



In nano-scale semiconductor devices the quantum effects become increasingly important 
and may dominate the transport phenomena. Whatever the system under study the tradi- 
tional semiclassical Boltzmann equation is no longer an adequate approximation. A rigorous 
quantum-mechanical approach is now necessary not only to study devices intrinsically based 
on quantum phenomena (e.g. resonant structures), but also to describe more conventional 
devices as nano scaled field effect transistors where quantum effects cannot be neglected 
anymore. However, the quantum transport theories are more difficult to implement within 
a general framework. Among the different quantum formalisms developed during the last 
decades,— the Non-Equilibrium Green Functions can be regarded as one of the most powerful 
due to its transparent physical route. It provides a general approach to describing quantum 
transport in the presence of scattering, including the contact and/or gate couplings. A very 
pedagogical review of the subject can be found in,— where the procedure of the self-consistent 
solution of the Green function and the Poisson equation is properly described. As stated by 
Christen and Biittiker,— the self-consistent treatment of the Coulomb potential is necessary 
to ensure gauge- invariant current- volt age (I-V) characteristics in non-linear transport. 

However, though the formalism is clear, the NEGF-calculation of transport quantities 
is often technically complicated. Except a few studies realizing the full two-dimensional 
calculation,— it is mostly developed for one- dimensional (ID) transport problems possibly 
coupled with the two- or three-dimensional description of electrostatics.^^ The resonant 
tunneling diode is a typical ID quantum structure which has been extensively studied both 
theoretically and experimentally. Although the basic structure of this device is simple, 
the development of a fully self-consistent NEGF-calculation likely to give reasonable I-V 
characteristics is still a subject of many efforts . 8 ! 9 ' 10 Alternatively, the Wigner function 
formalis m 11 ' 12 ' 13 or the self-consistent solution of the Schrodinger and Poisson equations^ 
have been also used to calculate transport properties of RTDs. 

The aim of this work is to present results of the fully self-consistent NEGF-calculation 
of typical electronic and transport quantities such as potential and electron density profile, 
transmission coefficient, I-V characteristics, and current noise spectral density in double 
barrier resonant tunneling structures. The Green functions are calculated exactly within 
the framework of the tight-binding model and self-consistently coupled with the Poisson 
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equation. The calculation procedure though standard is formulated in a simple form that 
can be readily applied to any ID nano-scale devices and can be easily extended to more 
complicated structures. Compared to other works in the literature, we obtain a good de- 
scription of resonant behavior for RTDs of different designs and at different temperatures. 
Although the NEGF-method has been used to calculate the noise in various structures , 15 ! 16 
to our best knowledge, the present work is the first attempt to use this formalism to show 
both sub-poissonian and super-poissonian shot noise in RTDs, which is now the subject of 
great attention, both theoretically ^ 18 i 19 i 20 i 21 i 22 and experimentally— i 24 i 25 i 26 i 27 

The paper is organized as follows. Section II describes the model, formulates the problem 
and the calculation method. In section III the numerical results of the potential and electron 
density profiles, the transmission coefficient, and the I-V characteristics are presented and 
discussed in comparison with those reported in the literature. Section IV is specially devoted 
to the noise calculation. Throughout the work, an attention is properly paid to the practical 
manner of calculation technique. 



II. METHOD AND FORMULATION 



Let us consider a double barrier resonant tunneling device connected to two infinite 
contacts, the left (L) and the right (R) (Fig. 1). We assume that each contact is characterized 
by an equilibrium Fermi function, fL(R) = l/[exp(E — ^l{r)) /ksT + 1], with a chemical 
potential, Hl(r)- The x-direction is chosen to be perpendicular to the barriers. The (y — z)- 
cross-section is assumed to be large so that the potential can be considered translation- 
invariant in this plane. The device is then described by the standard single-band effective 
mass Hamiltonian: 

which is coupled to two contacts. In Eq. (1), E c includes the bottom of conduction band 
and the conduction-band discontinuities, m\\(m±) is the longitudinal (transverse) electron 
effective mass, and U(x) is the total electrostatic potential. Within the framework of the 
single particle model, the potential U(x) can be determined from the Poisson equation: 

Ue d ^) = -[n d {x)-n e {x)\, (2) 
ax ax e 
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where e is the electron charge, eo is the dielectric constant of vacuum, e is the relative 
dielectric constant which may vary along the x-direction, according to the material, and n e 
and rid, respectively, are the electron and donor density, which are assumed to be independent 
of y- and z-coordinates. As usual, we consider the case of fully ionized donors. 

In order to solve the Hamiltonian (1) and Eq. (2) self-consistently with respect to the 
potential U and the electron density n e , taking into account the contact couplings, and 
further, to calculate the quantities of interest such as the electron density profile and I-V 
characteristics of the device, it is convenient to use the NEGF method. For the numerical 
procedure the system under study is spatially meshed and all quantities of interest are 
computed at the grid sites. For a grid spacing a along the x-direction, since the Hamiltonian 
(1) is already variable-separated, we can write it in the matrix form: 

[#]i,k,i',k' = — ~~ (2£b + U% + Ed + e(k)) S^i + ijA-iy] ^k,k'> (3) 

where z'(i') indicates grid sites, numbered from left (1) to right (N), t x = h 2 /2m\\a 2 is the 
coupling energy between adjacent sites along the x-direction (in the nearest neighbor tight- 
binding scheme), and e(k) = h 2 k 2 /2m± with k being the wave vector in the (y — z)-plane. 
The Hamiltonian matrix (3) is written in the {| | k >}-basis.— 

Once we have a matrix representation of the Hamiltonian operator it may seem straight- 
forward to get the retarded Green function (GF) by inverting the matrix, [(E — irj)I — H] -1 , 
where I is the unit matrix and r] is an infinitesimal positive real quantity, rj — ► + . In 
practical calculations, however, avoiding to work with infinite contacts, the GFs are calcu- 
lated only in the device domain and the contact couplings are introduced as the 'self-energy' 
matrices E r . Thus, for the device under study, taking into account the couplings with two 
contacts, L and R, the retarded GF can be determined as 

G r =[(E-i V )I-H-Tr L -T^-\ (4) 

where is the retarded self-energy matrix describing the coupling between the device 

and the left (right) contact. Within the framework of the ID nearest-neighbor tight binding 
model, these self-energy matrices can be calculated exactly: 

[ S L(i?)kky,k' =t 2 x [g L (R){E } \^]iy5^8 i)i i5i )1 (N) 1 (5) 

where gL(R)(E, k) is a solution of the equation: 

tl [9l(r)} 2 ~ <xl(r) (E, k)g L{R) + 1=0, (6) 
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with the convenient sign for the root. In this equation a^R){E,k) = E — e(k) + 2t x — 
([/(°) 4- Ec^)l(r), where, in the last term, E^, R -. is the bottom of conduction band at the 
L(R)-contact and U^) R s is defined by the applied bias. For instance, if the bias V is applied 
to the right contact, then = and U R ^ = V. 

The advanced GF and the advanced self-energies are Hermitian adjoints of the corre- 
sponding retarded matrices: G a = [G r } + and S a = [S r ] + . 

In NEGF method the lesser GF plays the central role, in term of which the measurable 
quantities are expressed. Noting that in the device under study the contact coupling is the 
only 'scattering' involved, the lesser GF can be expressed as^ 

G<(E,k) = G r (E,k)in(E,k) + T<(EMG a (E,k) 

= i[A L (E, k)f L (E) + A R (E, k)f R (E)]. (7) 

Here, we introduce the lesser tunneling rate defined as ^l(r) = ^l(r) ~ ^L(R)}fL(R)(E), 
and the spectral function Al(r) = G r [T l r L ^ — Y> a L ^ R ^\G a . As usual,— the equilibrium realized 
in a contact is assumed to be maintained even in the adjacent region of device. The Fermi 
function (/^ or fa) with the same chemical potential (/i^ or /zr) can be then applied in the 
corresponding region (close to L- or R-contact). 

Once the lesser GF (Eq. EJ) is known it is easy to calculate the electron density: [n e ]j = 
—iff dEdklG^ij. Assuming that the vertical effective mass m± is constant along the x- 
direction, the GF G K depends on k only through (E — s(k)), the sum over k- vector can be 
then changed into an integral that results in 

m ± k B T [^^[^ m\r^a-E, 

a=L,R 

where JF is the zero-order Fermi-Dirac integral, F^x) = ln(l + e x ). 

Thus, basically, the main body of the problem is to solve the Poisson equation (2) and 
to calculate the electron density (8) self-consistently. Once the self-consistent solution has 
been found, the current can be calculated as 

Ii = I = IoJ dET(E){T (^j^) - ^o(^f^)}. (9) 
Here Iq = (em±kBT)/(2iT 2 h 3 ) and the transmission probability matrix T(E) is defined as 

T(E) = T L (E)G r (E)T R (E)G a (E) , (10) 
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[dEj2 [AakmM^Y^y, i = 1, 2, N, (8) 



where = 1(1^,^ — ^1{R)) * s ^ e tunneling rate. 2 

Before presenting the numerical results we would like to mention that using this self- 
consistent treatment of the electron interaction, the NEGF approach used in this work 
becomes equivalent to the scattering Fisher-Lee theory.^ 2 - Furthermore, to solve the Pois- 
son equation we use the Newton-Raphson method with a Jacobian to be determined.— 
In this method it is unnecessary to find the Jacobian exactly and a good approxima- 
tion may be acceptable to reach the convergence. The Jacobian is actually estimated by 
fitting the electron density resulting from Eq. (8) with the electron density expression 
[n e ]i = N c J r i/ 2 ((fii — Ui)/ksT), where \ii is the fitting parameter, T\ji is the Fermi-Dirac 
integral of order 1/2, and N c is the effective conduction band density of states. Using this 
procedure, both Eqs. ([T]) and (0) are considered as a single non-linear one for U, which 
really reduces the time of calculation. 



III. NUMERICAL RESULTS 



The numerical calculations have been performed for typical 
Al x Gai_ x As / GaAs / Al x Ga>i- x As RTDs. Two barriers are symmetrical with the height of 
0.3 eV and the thickness of d = 3nm. The width of quantum well is w = 5nm. The barriers 
and the quantum well (scattering region) are undoped. The double-barrier structure is 
embedded between two low-doped buffer layers, each lOnm - thick. The whole system, in 
turn, is embedded between two outer access regions of GaAs (emitter and collector), doped 
at 10 18 cm~ 3 . The thickness of each access region is chosen as 30nm, which is believed to 
be large enough (see the inset of Fig. 4). The effective mass is assumed to be constant in 
the whole device and equal to 0.067m , where m is the free-electron mass. The relative 
dielectric constant of 11.5 is used. This system exhibits one resonant level associated with 
the conduction band minimum, approximately 0.14eU above the bottom of conduction 
band of GaAs, if the band bending is neglected. The temperature is generally chosen 
to be 300K. In particular, to check the device-size and the temperature effects, some 
calculations have been also performed for devices with quantum-well width of Anm, with 
barrier thickness of 2nm (Figs. 2 and 3), or at temperature of 77 K (Figs. 4 and 5). 

Here it should be noted that buffer layers with an appropriate thickness and dopant den- 
sity are generally used in numerical calculations as well as in experiments . 1 ?i 24 i 2 7 These layers 
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yield a reduction of the highest charge density at resonance in the quantum well and of the 
charge accumulation in the emitter region. Intrinsically such charge accumulations strongly 
affect the shape of the potential profile and therefore produce a current instabilit y 20 ! 28 How- 
ever in this work we limit our investigations to the structures with imperceptible instability. 

To model the Hamiltonian (3) the ID-grid of spacing a = 0.25nm is used. With the 
algorithm mentioned in the preceding section, typically, only 5 iterations are required to 
reach a self-consistent solution of Eqs. (2) and (8). Fig. 1(a) shows the potential profile cal- 
culated for the device with barriers 3nm- thick and the quantum-well 5nm-wide (hereafter, 
written for short as [3/5/3]-device) at various applied bias. At zero bias (solid line), the 
potential is lightly risen around the barriers that pulls the resonant level a bit upward. With 
an applied bias compared to the resonant one (dashed line) or higher (dot-dashed line), a 
large portion of the potential drop occurs at the right side of the barriers. Qualitatively, 
our results of potential profile are similar to those obtained in& for the [3/3/3]-device, using 
the stationary scattering matrix approach. A similar potential profile was also recently re- 
ported by Pinaud,— solving self-consistently the Schrodinger and the Poisson equations for 
the [5/5/5]-device. 

The electron density profiles, corresponding to the potentials in Fig. 1(a), are shown in 
Fig. 1(6). Clearly, for any bias, in the regions where the total potential is zero, the electron 
density is precisely given by the donor density. At zero bias the profile is certainly symmet- 
rical. A finite bias brings about an electron accumulation in the well. The accumulation 
reaches the highest level at the resonant voltage and decreases as the bias is continuously 
increased. At the same time, the profile becomes more and more asymmetrical. For the 
bias above resonance, a slight Friedel oscillation of the profile observed in the buffer layer 
adjacent to the left barrier is due to an electron accumulation at the left barrier, which 
occurs when the resonant level is much lower than the chemical potential in the left contact. 
Overall, these results are again in qualitative agreement with those presented in,— where 
however the data for the bias above resonance is not found. 

The most profound manifestation of all that are shown in Figs. 1 can be observed in 
Fig. 2(a), where the transmission coefficient T is plotted versus the energy E for the same 
device at different applied biases. For low bias (less than or about the resonant bias in I-V 
characteristics, V p ~ 0.311/, as can be seen in Fig. 5), the T(i?)-curve exhibits sharp peaks 
of almost unit on resonance and then falls off rapidly with energy on both sides. In the range 
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of energy under study we identify two resonant peaks in each T(E)-cmve separated by a 
distance weakly sensitive to bias and being equal to ~ 0.24eV. The bias, however, shifts 
the picture to the left making the first resonant peak to be cut off when the bias becomes 
higher than the resonant one (the case of 0.46V in the figure). Such a disappearance of the 
sharp peak describes the off-resonant state of the device. Certainly, both the width of peak 
as well as the peak-to-peak distance strongly depend on the device dimensions. In Fig. 2(6) 
we compare the T(£')-curves of three devices slightly different in the barrier thickness or 
the quantum- well width: [3/5/3] (solid line), [3/4/3] (dashed line), and [2/5/2] (dot-dashed 
line). In agreement with the data presented in,— Fig. 2(6) reasonably demonstrates that 
a decrease of either the quantum-well width or the barrier thickness makes the resonant 
peak wider and the peak-to-peak distance longer. In particular, two devices with the same 
quantum-well width of 5nm have the same position of the first resonant peaks, but the peak 
width is larger and the second peak locates at higher energy for the device with narrower 
barriers. 

In Fig. 3 we show the I-V characteristics for three devices corresponding with the trans- 
mission coefficients presented in Fig. 2(6). To analyze the data in Fig. 3 it is convenient to 
introduce two quantities: Ip is the current at the peak and 7 = Ipjly is the peak to valley 
ratio, where Iy is the current at the valley. Then, Fig. 3 demonstrates that a decrease of 
either the barrier thickness or the quantum-well width results in an increase of Ip, but a 
decrease of 7 ( 7 » 2.3/0.4, « 6/1.4, and « 8/3.4 for [3/5/3]-(closed circles), [3/4/ 3] -(black 
squares), and [2/5/2] -(black rhombus)-device, respectively). The latter statement is a con- 
sequence of the smearing of the resonant level that broadens the peak of T(E) as can be 
seen in Fig. 2(6). The width of resonance in energy is inversely proportional to the barrier 
thickness and/or the quantum-well width. That also explains, for example, why two devices 
with the same quantum-well width of 5nm, but with different barrier thicknesses, (3nm and 
2nm), exhibit not only strongly different Ip, but also different resonant biases, V p ~ 0.307V 
for [3/5/3]- (closed circles) and ~ 0.34V for the other, though the resonant energies are the 
same, i.e. « 0.14eV (see Fig. 2(6)). 

Concerning the temperature effect we show in Fig. 4 the I-V characteristics for the 
same [3/5/3]-device, but at different temperatures, 300K (closed circles) and 77fT(black 
squares). As the temperature is lowered from 300f^ to 77 K, the current Ip grows from « 
2.28 x 10 5 AcmT 2 to « 2.62 x 10 5 Acm~ 2 , and simultaneously, Iy falls from « 0.45 x l0 5 Acm~ 2 
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to ~ 0.33 x 10 5 Acm~ 2 , which finally produces a large change of the peak-to-valley ratio 7 
(from 7 w 5 to « 8). 

As an important note, in the inset of Fig. 4 we compare the I-V characteristics of 
two devices with the same quantum-well width and the same barrier thickness, but with 
different thicknesses of the GaAs access region, 30nm (closed circles) and 40nm (crosses). 
It is apparent that two curves are well coincident in the whole range of bias under study. This 
makes an argument to suggest that, as mentioned in the first section, the GaAs access region 
of 30nm thick can be considered as large enough and consistent with boundary conditions. 
Such an access region was used in calculations throughout this work. 

IV. NOISE CALCULATIONS 

In this section we demonstrate how the shot noise can be calculated using the NEGF- 
code developed in section II. Deviations of the noise from the full (Poissonian) noise value in 
RTD have been extensively investigated in a great number of works, both theoretica l 19l2 ° l22 
and experimental . 23 ! 24 ^ 25 ! 26 ^ 27 Mathematically, the measure of these deviations is the Fano 
factor F defined as the ratio of the actual noise spectral density to the full shot noise 
value 2eJ, where I is the average current. Physically, it is widely accepted that the Pauli 
exclusion and the charge interaction are the two correlations, which cause observed shot 
noise deviations. While the Pauli exclusion always causes a suppression of noise, the charge 
correlation may suppress or enhance the noise, depending on the conduction regime. In 
RTD, is was experimentally found that the noise is partially suppressed (sub-Poissonian 
noise, F < 1) at low bias voltages (pre-resonance) and becomes very large (super- Poissonian 
noise, F > 1) in the negative differential conductance (NDC) region. From theoretical point 
of view, typically, there are two approaches based on different descriptions of tunneling 
process. 

Treating the tunneling through RTD as a quantum coherent process, a general expression 
of noise has been derived . 17 i 20 In the limit of zero frequency and for two-terminal structures 
it has the form: 



S 




(11) 
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where the first term, proportional to f a (E)(l — f a (E)), describes the thermal noise, and 
the second term, proportional to T(E)[1 — T(E)], is nothing but the partition noise. This 
expression leads to an important consequence that the Fano factor should never be less 
than 1/2. With T(E) determined in section II, the noise (11) can be straightforwardly 
calculated. In the inset of Fig. 5 the Fano factor F calculated in this way for the [3/5/3]-device 
is presented in the range of low (pre- resonant) bias V. It is clear that the obtained F(V)- 
dependence describes quite well existent experimental data , 19 ' 24 including a noise suppression 
in the positive differential conductance (PDC) region close to the resonance. At higher 
bias, in the NDC region, as critically discussed in l^jj, the charge interaction becomes to 
dramatically affect transport properties of the quantum well, and the expression (11) can 
no longer be used. 

On the other hand, treating the tunneling as sequential process in the spirit of master 
equation, Iannaconne et a/ . 18 ' 19 arrived at the following expression formula for Fano factor: 

F = 1 - 2r g T r /{r g + r r ) 2 , (12) 

where the generation and recombination rates through barriers are given by 

r" 1 = -(dI { L +) /dN) \ N ., r- 1 = (dI { - ] /dN) \ N ., (13) 

N is the number of electrons in the well at given bias point and N* is the steady-state value 
of N, defined from the charge conservation l£ (iV*) = I R ~\N*). The quantities 1^ and 
I R are here the flux of electrons injected into the well through the left barrier and the flux 
of electrons leaving the well through the right barrier, respectively. It was reported thalj 1 ^ 
the Fano factor calculated from (12) is in good agreement with the experimental data for the 
[12.4/6.2/14.1]-sample in a large range of bias, including both the PDC pre-resonant region 
with a sub-poissonian noise and the NDC region with a super-poissonian noise. To realize 
a calculation of noise (12), using the NEGF-algorithm developed in section II, we note that 
the fluxes 1^ and I R ~^ can be actually expressed as 

r(+) _ f AT? F. f T?\\ A . I T?\ I A - I FM %( ^ L — — ^ 



k R T 



4 +) = J dE T L (E)[A L (E) + A r {E)]Tq 

4 +) = JdE T R (E)[A L (E)F (^j-f) + A R (E)M^j-^)]., (14) 

where Al(r) and ^{l{r) are already defined in eq.(7) and eq.(10), respectively. So, actually, 
we can also calculate the Fano factor F (12), using the NEGF-code developed. In Fig. 5 we 
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show F calculated in this way as a function of bias for the [3/5/3] -device at two temper- 
atures, 77 K (closed circles) and 300.K" (black squares). Overall, the obtained F(l/)-curves 
qualitatively resemble experimental dat a 19 ' 24 with a minimum at the resonant bias and a 
super-poissonian peak in the NDC region. For a given device, with decreasing temperature, 
as can be seen in Fig. 5, the minimum becomes deeper (still greater than 1/2) and the super- 
poissoinian peak becomes higher. The dramatic increase of noise observed in experimen t 19 ^ 24 
is caused by the charge interaction enhanced by the particular shape of the density of state 
in the well. In deriving the expression (12) such an interaction has been self-consistently in- 
cluded. The role of self- consistent Coulomb interaction in creating a super-Poissonian noise 
is also recently demonstrated by Oriols et al.. 22 

Thus, using the NEGF-algorithm developed in section II we are able to calculate the cur- 
rent noise in both tunneling regimes, coherent and sequential. It should be here mentioned 
that in literature there are different opinions about how a super-Poissonian noise can be 
produced and its relation to the NDC. In the theory^ 9 - the super-poissonian noise is closely 
associated with the the NDC, when the time r g (12) is negative. Blanter and Buttiker em- 
phasized the role of current instability and stressed that the bias range of super-poissonian 
noise is generally different with the NDC range. 20 Comparing the I-V curves and the noises 
measured in a super-lattice diode and in a RTD, Song et alM concluded that not NDC, but 
charge accumulation in the well responds for the super-poissonian noise observed in RTD. 
Safanov et a/.—, measuring the noise in resonant tunneling via interacting localized states, 
observed a super-Poissonian noise in the range of bias, where there is no NDC. Authors have 
also pointed out that the effect on noise of the Pauli exclusion principle and of the Coulomb 
interaction are similar in most mesoscopic systems. For Coulomb blockade metallic quantum 
dot structures, studie d 30 ' 31 support the idea that not NDC, but charge accumulation in dots, 
responds for super-poissonian noise. Remarkably, Aleshkin et alM recently claimed that 
while in the sequential tunneling double barrier model the Fano factor F is still limited by 
the lowest value of 1/2, in the coherent model F may drop below this value. Then, authors 
also suggest that the noise suppression below the value 1/2 can be used to identify a coherent 
transport. Actually, to solve all these, very fundamental, contradictions a more systematic 
analysis of noise, taking adequately into account the charge interaction is requested. To this 
end, it is perhaps most convenient to use the NEGF, well-known as the quantum transport 
approach capable even to deal with far from equilibrium systems. 
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V. CONCLUSION 



The NEGF approach has been formulated and implemented in a fully self-consistent cal- 
culating procedure that can be readily applied to any ID nanoscale structures and extended 
to more complicated devices, e.g. nanoscale field effect transistors. The Poisson equation 
solver routine has been improved to considerably speed up calculations. It yields physically 
reasonable results and allows us to work in a large range of temperature. 

Numerical calculations have been performed for RTDs of different designs and at different 
temperatures. The potential and electron density profiles obtained for various applied bias, 
below and above resonance, rationally describe the resonant behavior in the device. The 
transmission coefficient and the corresponding I-V characteristics seem to be sensitive not 
only to the quantum-well width, but also to the barrier thickness. Besides the well-known 
fact that the quantum-well width defines the resonant level, our calculations show that 
a decrease of either the quantum-well width or the barrier thickness makes the peak of 
transmission coefficient wider and the peak-to-peak distance longer. Correspondingly, the 
peak-to-valley ratio in I-V curves decreases, though the value of the current at the peak is 
considerably risen. The reason merely lies in a broadening of the resonant level, caused by 
narrowing the quantum- well width and/or the barrier thickness. Additionally, by lowering 
the temperature from 3007^ to UK we observed not only a raise of the peak in I-V curve, but 
also a reduction of the valley current that leads to a significant increase in the peak-to- valley 
ratio. 

It was also shown that the NEGF-algorithm developed can be straightforwardly used to 
calculate the current noise in both coherent and sequential tunneling models. In qualitative 
agreement with experiments, obtained results highlight the role of charge interaction which 
causes a fluctuation of density of states in the well and therefore a noise enhancement in the 
NDC region. We believe that in this way, i.e. using the NEGF, a fully quantum transport 
approach for non-equilibrium problems, one can calculate the noise in a variety of interacting 
mesoscopic systems and therefore better understand the nature of noise deviations as well 
as of conduction mechanisms in these structures. 
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FIG. 1: Potential profile (a) for the [3/5/3] -device and corresponding charge density (b) at various 
applied bias: OV (solid line), 0.31F (close to resonance, dashed line), and 0.46F (above resonance, 
dot-dashed line), Temperature: 300K. 

FIG. 2: (a) Transmission coefficient in logarithmic scale for the [3/5/3] -device at 300K at the 
applies bias: 0.31V^ (dashed line), and 0.461/ (dot-dashed line), (b) Comparison the transmission 
peak of different structures: [3/5/3] (solid line), [3/4/3] (dashed line), and [2/5/2] (dot-dashed 
line) at V = 0.00V. 

FIG. 3: I-V characteristics for the same devices as in Fig. 4: [3/5/3] -(closed circles), [3/4/3] -(black 
squares), and [2/5/2]- (black rhombus); Temperature: 300K. 

FIG. 4: The temperature effect: 300i^(closed circles) and 77K (black squares) for the same 
[3/5/3]-device. Inset: the I-V characteristics for two devices with the same quantum-well width 
and the same barrier thickness, [3/5/3], but with different thicknesses of the GaAs access region, 
30nm(») and 40nm(x), are compared. Temperature: 300K. 

FIG. 5: The Fano factor calculated from Eq.(13) is plotted as a function of the applied bias for the 
[3/5/3] -device at 300K (black square) and 77K (closed circles). Inset: the Fano factor calculated 
from Biittiker's formula (11) for the same device at 77K. 
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FIGURE 1 (V. Nam DO) 
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FIGURE 2 (V. Nam DO) 
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FIGURE 3 (V. Nam DO) 
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FIGURE 4 (V. Nam DO) 
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FIGURE 5 (V. Nam DO) 
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